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ABSTRACT _ 

This paper discussed how wind farm disturbance, especially wind speed 
variability, can affect the performance of the power system. Also, it 
discussed how blades angels of wind turbine can be controlled to increase the 
energy efficacy of the power system. It showed that the optimized pitch angel 
controller using harmony search algorithm could enhance blades angels’ 
adjustments performance. First part, the paper explained the advantages of 
doubly fed induction generator in wind turbine system. Paper also 
enumerated the most probable ways to tame the wind speed variability 
challenge focusing on the pitch angle controller technique. After that, paper 
compared the system parameter's result before optimization with values after 
optimized pitch angle controller gain (Kp). This comparison would be held in 
three cases, case of variable wind speed with normal operation condition, 
constant wind speed with line to ground fault condition and variable wind 
speed with line to ground fault condition. Finally, It demonstrated the 
MATLAB/SIMULINK model used illustrating results appeared 
and conclusion. 
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1. INTRODUCTION 

Increased growth of energy demand coupled with fossil fuel resources depletion and environment 
concerns has urge society to search for an alternate renewable energy sources. Out of available renewables, 
Wind energy has proved itself as an attractive solution for wide geography for many reasons [1]. It's 
considered as agreen energy source, it needs low running cost and it allows more than 90% of used land to be 
used for cultivation. 

Fixed speed wind turbine (FSWT) and variable speed wind turbine (VSWT) are the most popular 
types of wind turbines as they are able to extract optimum power which improving overall efficiency. 

One of the best variable speed wind turbine is the doubly fed induction generator (DFIG). DFIG is 
currently one of the most popular wind energy extraction system Figure 1 [2-4]. DFIG solution has proved 
itself as one of the most practical wind energy extraction solution. With reduced Converters rating 
approximately (20-30) % of the total generator power. This means that the losses in the power electronic 
equipment can be reduced in comparison to FSWT or VSWT. Speed variation of (30%) around synchronous 
speed can be obtained which means more energy captured from the wind as shown in Figure 2. Also it can 
transmit power for the grid through both stator and rotor. Moreover, it offers reduced mechanical stress on 
the wind turbine [5]. The DFIG allows decoupled control of active & reactive power [3, 6, 7]. 
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Figure 1. Double fed induction generator (DFIG) 
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Figure 2. Comparison between different turbines for Produced Energy at different wind speeds 


2. WIND SPEED VARIABILITY 

Stochastic nature of wind speed comes from its dependency on weather and climatic conditions. 
This comprises a crucial Challenge for Power Extraction operation [2, 8]. Wind Speed variability is 
considered the main challenge for many power system issues including Power quality, Power factor and 
voltage regulation. 

Figure 3 shows a particular site in Ireland, A (36-hour) ahead wind speed forecast in hourly steps for 
three distinct days is done as shown in Figure 3, making a comparison between the actual and forecasted 
values. In general, wind speed predictions are seen to be reasonable, with little fall-off in performance for 
longer time horizons. Figure 3b illustrates fairly unusual day, when looking (lOhours) ahead the actual wind 
resource is significantly greater than predicted [8]. 





Figure 3. 36-hour ahead wind speed prediction (a, b and c) 
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Wind speed variability issue and its effect on power system quality have motivated reaserchers to 
obtain asolution to this problem. There are many techniques have been studied to overcome this problem. 
Amoung these solutions, rotor and stator side converters control are produced [9-10]. Also, pitch angle 
control technique has proved itself as one of the most effective wind speed control techniques enhancing 
wind turbine performance. It has been studied using many optimization techniques such as Fuzzy logic 
control [11]. Herein, it will be studied based on harmony search algorithm. 


2.1. Pitch angle controller 

Pitch angle controller enables operator to adjust all blades at the same time or each one separately. 
By pitching the blades towards wind maximizing energy captured or pitching blades away discarding the 
excess power, the blades regulate the power captured by rotor keeping power and speed at their rated values. 
Also, to prohibit any mechanical failure due to blades over speed [4]. Pitch angle control is limited to about 
(107s) during normal operation and (207s) for emergencies [3]. The extracted wind energy is affected by the 
rotor power coefficient (Cp) as in [1], [2]. Determination of that coefficient is dependent on the aerodynamic 
parameters for each turbine. Performance curves which provide information about efficiency of the turbine 
for different wind speeds and various blades angles can also illustrate the condition at which maximum 
power coefficient (Cp) could be obtained. This is shown in Figure 4 [12]. To calculate it the following 
equations are used [13]. 

P t =\pv 3 C p ( 1 ) 




Pwind turbine 

Pair 
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(3) 



1 

a+c 6 /? 


—) 
P 3 +1 J 


(4) 


Where Pt is the mechanical output power of the turbine (w). p is the air density (approximately 
1.225 kg/m 3 ). A is the swept area of rotor (m 2 ). v is the wind speed (m/s). Cl, C2,..., C7 are the turbine 
Aerodynamic, p is the pitch angle. X is the turbine tip speed ratio. Form (3-4) Performance curve can be 
drawn for different pitch angle values.as shown in Figure 4. These curves illustrate the maximum energy 
extraction points and show that the power coefficient is inversely proportional to the pitch angle value. 



Figure 4. Power coefficient vs. tip speed ratio for various values of pitch angle 


Figure 5 shows a model of pitch angle controller. A reference rate of change of pitch is generated 
using PI controller. This reference is limited. A first order system gives the dynamic behavior of speed 
control of pitch variation. By integrating the variation of the angle, pitch angle is obtained [3]. 

In this paper, pitch angle control based on harmony search optimization technique will be 
considered. System in Figure 6 will be studied at three cases. The first case will be the system at normal 
condition, in case of variable wind speed. Second one will be the system at fault condition, line to ground 
fault, in case of constant wind speed. Finally, the system at fault condition, line to ground fault, in case of 
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variable wind speed. In each case optimization of blade angle control will be carried out by the harmony 
search algorithm (HAS) techinique. 
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Figure 5. Pitch angle controller configuration 


3. DFIG WIND FARM MODEL IMPLEMENTATION AND SIMULATION RESULTS [14] 

Figure 6 shows a wind farm represented by one (1.5 MW) wind turbine connected to a 120 kV, 60 
HZ frequency power system through (2500MVA) Transmission line. A step down transformer (120/25) KV 
of 47 MVA power is connected to a 25 Kilo- meters transmission line, where the studied fault is subjected in 
the middle of it. Then to a step down transformer ( 251.515 ) KV of 1.5 MVA power through appoint of 
common coupling with the wind turbine and an arbitrary load of 20 watt. A doubly-fed induction generator, 
DFIG, which is consisting of a wound rotor induction machine. Its stator is directly connected to the power 
grid while rotor is connected to power grid through AC/DC/AC IGBT-based PWM converters. A DC Bus 
link of 1200V and Capacitor of size (10000e-6) farad. With starting Pitch angle gain value of (Kp=500) and 
maximum pitch angle of 45 degree [14]. 



Figure 6. DFIG wind farm model represented by MATLAB/SIMULINK 


3.1. Harmony search algorithm [15,16] 

It’s a music-based meta-heuristic optimization algorithm. It was inspired by the improvisation of 
Jazz musicians as shown in Figure 7. IPs firstly proposed by Zong Woo Geem in 2001. Figure 7 shows the 
analogy between the music improvisation and the engineering optimization process. 
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Figure 7. Analogy between the Music improvisation and engineering optimization process 


In HSA, each decision variable generates a value to find the best harmony, global optimum, all 
together as shown in Figure 8. HS algorithm does not require initial values for the decision variables. Also, 
HS uses stochastic random search instead of gradient search. This random search is based on the harmony 
memory considering rate and pitch adjusting rate so that derivative information is not necessary [15-16]. 


Step 1 Step 2 



Figure 8. Harmony search algorithm optimization process 


3.2. Objective function [17]: 

The objective function aims to minimize the error of voltage and active power. This means better 
system performance, faster response and longer connection time between wind turbine and grid. Many 
performance criterias are used to reach optimum solution. For example; integrated absolute error (IAE), 
integrated square error (ISE), Integrated time-weighted-squared-error (ITSE) and-integral time-weighted- 
absolute error (ITAE). All of these methods are represented in the next [5-8] respectively. 


IAE = J 0 f | AV e |dt (5) 

ISE = f*(AV e ) 2 dt (6) 

ITSE = f* t(AV e ) 2 dt (7) 

ITAE = j* t\AV e \dt (8) 


The used objective function is ISE in all calculations. It's characterized by small maximum 
overshoot (Mp) and long settling time (Ts) as it weights all errors equally, time independent. Although ITSE 
method can overcome his demerit, but it has complex analytic formula and longer time consumption. 
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3.3. Simulation results and analysis 

To find out the effect of blade angle control on power system stability of grid connected wind 
turbine. The System shown in Figure 6 will be studied at three different cases. Firstly, it will be studied at its 
normal operation condition, with no fault occurrence, in case of variable wind speed. Secondly, it will be 
studied at fault condition, line to ground fault, in the middle of the transmission line at constant wind speed. 
Thirdly, it will be studied at its faulted condition, same as previous case study, but with variable wind speed. 
A comparison will be done in all previous case studies with the optimized gain (Kp) of the pitch angle 
controller. Pitch angle controller based on HAS will be used. The Line to ground (L-G) fault in the middle of 
transmission line is supposed to be between (20 to 25) seconds. 

3.4. First case, normal operation and variable wind speed (before & after) optimization 
3.4.1. Wind speed curve:- 

Figure 9 shows wind speed curve illustrating its variability nature. Due to wind speed variance 
especially at (sec 27 till sec 44.11) a reduced energy is captured by blades. This may lead to under speed 
condition. This wind speed will be applied in model of Figure 6 in the next sections studying its effect on 
voltage and power curves performances with and without optimized blade angle controller gain “Kp” as 
shown in Table 1. 


i 

1 - 


Figure 9. Wind speed curve versus time 


3.4.1.1 Voltage at bus B575 (P.U) 

Figure 10 shows voltage curves at bus (B575) before and after optimization at their rated value 
(1 P.U). There are no changes noticed in voltage curve as there is no fault occurred. The wind speed varience 
will affect delivered power to the generator without affecting voltage of the bus at common coupling 
point (BCC). 
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Figure 10. Voltage (before & after) optimization at wind farm Bus (B575) 
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3.4.1.2 active power (P.U) 

Figure 11 shows active power curves (before & after) optimization. Figure 11 shows preference of 
power curve after optimization over its values without optimization. It shows that the optimized gain of pitch 
angle has increased energy captured enhancing energy efficiency. Also, it prevented power curve from falling 
at time between (sec 16.75 to sec 20.22), (sec 21.5 to sec 23) and (sec 44.5 to sec 45.2) at which the power, 
before optimization, has fallen to values of (0.2 P.U), (0.8 P.U) and (0.73 P.U) respectively due to reduction 
in wind speed. These optimized values of blades angles helped the wind turbine to achieve stable and reliable 
operation. Also, it enhaned its ability to stay connected to the grid. 





Figure 11. Active power curves (before & after) optimization 

3.4.1.3 Reactive power (P.U) 

Figure 12 shows reactive power curves (before & after) optimization. It shows a slightly difference 
between two curves. At time between (sec 17 to sec 20.3) reactive power, before optimization, has been 
injected to the system to compensate reduction in delivered power. The reactive power curve, after 
optimization, shows more stable operation as optimized pitch angle controller. At time between (sec 29.5 to 
sec 44) both curves shows rising up to inject the system with reactive power as the wind speed sharply 
reduced at this time. 
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Figure 12. Reactive power curves (before & after) optimization 


3.4.1.4 Pitch angle (deg) 

Figure 13 and Table 1 show pitch angle curves (before & after) optimization conditions. Figure 13 
illustrates slightly changes between two curves especially at time between (sec 17 to sec 23.3). It shows a 
different behavior of optimized blades angles in order to keep energy captured close to its rated value (0.9 
P.U) despite of reduction in wind speed condition. 
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Figure 13. Pitch angle curves (before & after) optimization 


Table 1. Pitch Ang le Controller Gain “Kp” values with and wi thout Optimization 


Kp without Optimization 

Kp with Optimization 

Kp = 500 

Kp = 97.0593 


Second case, fault condition, line to ground fault at the middle of transmission line, and constant 
wind speed (before & after) optimization. The fault occurrence time will be from (sec 20 to sec 25). The 
wind speed will be kept constant at (15 m/sec). 


3.4.2. Wind speed curve 

Figure 14 shows wind speed curve which is set to constant value of (15 m/sec). This case shall help 
for better understanding for the effect of optimized pitch angle controller gain (Kp) on both active and 
reactive power curves performance at fault condition. This wind speed will be applied in model of Figure 6. 
In the next sections studying its effect on voltage and power curves performances with and without 
optimization of blade angle controller gain “Kp” as shown in Table 2. 



Figure 14. Wind speed curve versus time 


3.4.2.I. Voltage at Bus B575 (P.U) 

Figure 15 shows voltage curves at bus (B575) before and after optimization. Figure illustrates 
aslightly improvement in voltage curve especially during fault. The voltage curve in case of optimized pitch 
angle controller has fallen to (0.9 P.U) while it has fallen to (0.8 P.U) in case of without optimization. 
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However, voltage curve after optimization is subjected to higher degree of oscillation than other curve. These 
oscillation can be decreased by injecting an improved controller technique into the system. 



Figure 15. Voltage (before & after) optimization at wind farm Bus (B575) 


3.4.2.2. Active power (P.U) 

Figure 16 clarifies the effect of optimized pitch angle on keeping the wind turbine connected 
through accidental disturbance such as line to ground fault. Figure 16 shows both active power curves (before 
& after) optimization rising almost to their rated value of (1 P.U) after overcoming the inertia forces and 
mechanical losses. Then, power curve before optimization has collapsed at time (8.3 sec) to zero value. This 
happened because of the wind turbine protection threshold which cut-off the turbine as it reached over 
current condition. On the other hand, active power curve after optimization oscillates above the rated value 
(between 1.5/0.5 P.U) values at time between (sec 8 to sec 14). Then, the curve settled down till the time of 
fault occurrence. It's shown that the power curve has raised at fault time (sec 20 to sec 25) helping wind 
turbine to stay connected and ride through fault. The power curve of optimized pitch angle is subjected to 
higher degree of oscillation which it can be reduced by injecting an improved controller technique into the 
system. 


1 

I 






Figure 16. Active power curves (before & after) optimization 


3.4.2.3. Reactive Power (P.U) 

Figure 17 shows reactive power curves (before & after) optimization. It shows reactive power curve 
before optimization at zero value till the turbine is tripped out from grid at time (8.3 sec) due to over speed 
condition. Also, the figure shows the after optimization curve of reactive power increased at fault time 
between (sec 20 to sec 25) to value of (0.9 P.U). Reactive power is injected at fault time to compensate 
current absorbed at faulted Point. Figure 17 declared the effect of optimized gain (Kp) of pitch angle in 
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keeping the wind turbine connected to the grid and riding through fault occurred. However, it shows high 
oscillation in reactive power curve, after optimization case, which is can be diminished by injecting an 
improved controller technique into the system. 



Figure 17. Reactive power curves (before & after) optimization 


3.4.2.4. Pitch angle (deg) 

Figure 18 and Table 2 show pitch angle curves (before & after) optimization. More controllable 
variation in blades angles has been indicated in pitch angle curve after optimization (Kp=13.8624). This shall 
enhance its stability. Also, it shows that lower values of pitch angles allowed more energy captured by wind 
turbine improving its performance. Optimized pitch angle gain (Kp) prohibited the power curve from 
collapsing to zero at the fault time helping the wind turbine to stay ON for longer time. This will help wind 
farm operator to ride through system faults and to fulfill the grid code requirements. 



Figure 18. Pitch angle curves (before & after) optimization 


Table 2. Pitch Ang le Controller Gain “Kp” Values with and wi thout Optimization 


Kp without Optimization 

Kp with Optimization 

Kp = 500 

Kp = 13.8624 
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Third case, fault condition, line to ground fault at the middle of transmission line, and variable wind 
speed (before & after) optimization. The fault occurrence time will be from (sec 20 to sec 25). 

3.4.3. Wind speed curve 

Figure 19 shows wind speed curve illustrating its variability nature. Due to wind speed variance 
especially at (sec 27) which is concurrent with pitch blades angle come back to last position due to its inertia. 
This leads to (under speed) condition which will lead to cut off the turbine by the wind turbine protection 
threshold. This wind speed will be applied in model of Figure 6 in the next sections studying its effect on 
voltage and power curves performances with and without optimization of blade angle controller gain “Kp” as 
shown in Table 3. 



•MM («#•) 


Figure 19. Wind speed curve versus time 


3.4.3.I. Voltage at Bus B575 (P.U) 

Figure 20 shows visional comparison between voltage performance before and after pitch angle 
controller gain (Kp) optimization. It shows both curves values at its rated values (1 P.U) till the fault instant 
(20 sec) when they reduced to nearly (0.87 P.U and 0.88 P.U) of their values in case of without and with 
optimization respectively.This (15%) voltage drop is because of passing the fault current at the faulted point. 
Voltage drop is slightly improved after optimization as the active power transmitted is prevented from 
collapsing to zero. Curve oscillation is clearly noticed in case of after optimization. However, the curve 
oscillation can be reduced by injecting an improved controller technique into the system. 





Figure 20. Voltage (before & after) optimization at wind farm Bus (B575) 


3.4.3.2. Active Power (P.U) 

Figure 21 shows visional comparison between active power performance before and after pitch 
angle controller gain (Kp) optimization. After optimization; active power delivered is prohibited from falling 
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to zero value at fault time thanks to optimized pitch angle controller. Then the power, after optimization, is 
oscillating at fault clearance (sec 25). At time (sec 30) it settled at its rated value (1 P.U). Figure shows the 
active power before optimization is tripped off at (37.13 sec) while it lasts to longer time of (42.09 sec) in 
case of after optimization. Therefore, the optimized gain (Kp) has helped the wind turbine to be connected 
for longer time when the wind speed has reduced to low value close to the (under speed limit) threshold. The 
power curve in case of after optimization suffers from high oscillation which can be reduced by injecting an 
improved controller technique into the system. 



Figure 21. Active power curves (before & after) optimization 


3.4.3.3. Reactive power (P.U) 

Figure 22 shows reactive power curves (before & after) optimization. It shows that both curves are 
at their zero value till the line to ground fault occurred at time (20sec). At that time reactive power is injected 
through the system to retrieve absorbed fault current at the faulted point till the fault clearance time (25 sec). 
Figure 22 illustrates slightly improvement of curve and longer time of connection, but it shows a high degree 
of oscillation. These oscillations can be reduced by injecting an improved controller technique into the 
system. 



Figure 22. Active power curves (before & after) optimization 


3.4.3.4. Pitch angle (deg) 

Figure 23 and Table 3 show visional comparison between Pitch Angle Performance before and after 
optimization.The curve became lower, more stable as the variance of wind speed became more controllable. 
Thus helped the wind turbine to stay connected for longer time, from (37.13 sec to 42.09 sec), before the 
wind protection threshold trip it as it reached the Under Speed limit due to wind speed variance. 
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Figure 23. Pitch angle curves (before & after) optimization 


Table 3. Pitch Ang le Controller Gain “Kp” Values with and wi thout Optimization 


Kp without Optimization 

Kp with Optimization 

Kp = 500 

Kp = 4.6171 


4. CONCLUSION 

In this paper the Preference of DFIG over other machines have been discussed. Also an attempt is 
made to use the MATLAB/SIMULINK to investigate the effects of wind speed variability nature on Wind 
Turbine performance. Hence, it discussed how to control wind speed variability with the help of optimized 
pitch angle controller based on harmony search algorithm. As observed, simulation results indicate the 
influence of blades angels’ adjustment on the active power curve. It prevented it from falling to zero value 
helping the wind farm operator to ride through faults without being disconnected from the grid. Also it 
observed that pitch angle curve with optimized gain has become more smooth than its previous value without 
optimization .This would help the system to stay grid connected for longer period and ride through fault 
occurred in the middle of the transmission line even though the sharp variance of wind speed curve. As it 
reached the Under Speed limit determined by the wind turbine protection threshold, system has been tripped 
out from grid. The results of second and third cases illustrate that the system after optimization is subjected to 
high curve oscillation. These oscillations can be treated using any reactive power compensation technique or 
by injecting an improved controller technique into the system. 
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